Lambda-Perturbations of Keplerian Orbits 
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To estimate the influence of the dark energy on the planetary orbits, we solve the general rel- 
ativistic equations of motion of a test particle in the field of a point-like mass embedded in the 
cosmological background formed by the Lambda-term. It is found that under certain relations be- 
tween three crucial parameters of the problem — the initial radius of the orbit, Schwarzschild and 
de Sitter radii — a secular perturbation caused by the Lambda-term becomes significant, i.e. can 
reach the rate of the standard Hubble flow. 

PACS numbers: 04.25.Nx, 95.10.Ce, 04.80.Cc, 95.36.+X 
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The question if the planetary orbits and other local ce- 
lestial objects are subject to the cosmological influences 
(in particular, if they feel the universal Hubble expan- 
sion) was put forward by McVittie as early as 1933 [l|; 
and this problem attracted attention of a number of other 
researchers during the few subsequent decades I2 H12I . A 
quite comprehensive review was given by Bonnor [13j. 

The most frequent conclusion of these studies was that 
the effect of cosmological expansion at the planetary 
scales should be very small or absent at all. However, 
the particular estimates given by different authors sub- 
stantially disagree with each other. Moreover, the most 
of these estimates (excluding the recent ones) are not 
applicable to the case when the cosmological background 
is formed by the "dark energy" (i.e. the A-term in Ein- 
stein equations), just because it is distributed perfectly 
uniform and insensitive to the local gravitational pertur- 
bations. 

For example, the most well-known argument against 
the local Hubble expansion is the so-called Einstein- 
Straus theorem [2\: Let us consider a uniform back- 
ground cosmological matter distribution and, next, cut 
out a spherical cavity and concentrate all its mass in the 
central point. Then, a solution of the General Relativity 
equations will be given by the purely static Schwarzschild 
metric inside the cavity, and it will transform to the time- 
dependent Friedmann-Robertson- Walker metric outside 
the cavity. In other words, there is no Hubble expansion 
in the local empty neighborhood of the point-like mas- 
sive body, but such an expansion appears in the regions 
of space filled with the cosmological background matter. 
Unfortunately, despite an apparent generality of this re- 
sult, it is evidently inapplicable to the A-dominated cos- 
mology, because it is meaningless to consider an empty 
cavity in the vacuum energy distribution. Similarly, it 
can be shown that such arguments against the local Hub- 
ble expansion as the "virial criterion" of gravitational 
binding and Einstein-Infeld-Hoffmann surface integral 
method [J] also do not work when the cosmological back- 
ground is formed by the perfectly uniform A-term. 

Therefore, the most straightforward and self-consistent 
way to estimate how much can the dark energy affect the 



planetary dynamics is just to solve the two-body prob- 
lem in the pure A-background. In the simplest case of 
a test particle of infinitely small mass moving around 
the point-like mass M, this can be done using the well- 
known solution of General Relativity equations obtained 
long time ago by Kottler LJ] (in the modern literature, 
it is often called Schwarzschild-de Sitter solution): 



ds 2 = - 1- 



2GM Ar 



2GM Ar 



c 2 dt' 2 



dr" 



r'\d9 2 - 



(1) 
dtp 2 ) ; 



for a more general review, see also [15[. Here, G is the 
gravitational constant, c is the speed of light, and primes 
designate the original Kottler's "static" coordinates. 

Since metric (fTJ) was derived well before the birth of 
the modern cosmology, it suffers from a lack of the ad- 
equate cosmological asymptotics at infinity; namely, it 
does not reproduce the standard Hubble flow. Unfortu- 
nately, this fact was ignored in a number of recent stud- 
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24j . To resolve the above problem, it is necessary to 



perform transformation to the commonly-used cosmolog- 
ical Robertson- Walker coordinates (represented below by 
the variables without primes): 
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as outlined in our earlier work [l( 
metric will take the form: 
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where, for conciseness, we put c = 1, and dot denotes a 
derivative with respect to the proper time of the moving 
(4b) particle. 

Next, if the coordinate system is oriented so that the 
particle moves in the equatorial plane, 6 = ir/2 = const, 
these equations are reduced to the following set: 
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In the above formulas, r g = 2GM/c 2 is Schwarzschild 
radius, r\ = \/3/A is de Sitter radius, and ao is the 
scale factor of the Universe. 

Taking ao = 1 at t ~ and keeping only the lowest- 
order terms of r g and 1/^a, metric (|4a j) -(|4d" ]) can be 
rewritten in a more compact form: 
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Motion of a test particle in this metric is described 
by the geodesic equations which can be derived by the 
standard way: 
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In the case of a realistic planetary system, the prob- 
lem under consideration involves three characteristic 
scales, which differ from each other by many orders of 
magnitude — Schwarzschild radius r g (e.g., ^10~ 2 m for 
the Earth as a central body), a typical initial radius of 
the planetary orbit Ro (e.g., ^10 9 m for the Moon mov- 
ing around the Earth), and de Sitter radius ta (~10 m, 
which clepends on the amount of dark energy in the Uni- 
verse) [25j |. 

Despite of the availability of the small ratios, the an- 
alytical treatment of set of equations (j7a|) - (j7c|) is very 
hard, because the corresponding expansions in terms of 
the small parameters do not converge. The numerical 
integration is also challenging, because the standard ac- 
curacy of representation of the floating-point numbers 
in a computer is usually insufficient to cover the above- 
mentioned range of numbers. As a result, it is neces- 
sary to use a special software for emulation of the high- 
accuracy arithmetic. 

A detailed numerical analysis of these equations will be 
published elsewhere; while here we restrict our consider- 
ation by a few toy models, where difference between the 
characteristic scales is not so much as in reality. This will 
help us to reveal the most important features of the re- 
sulting motion. From here on, it is convenient to take the 
initial radius of the orbit Ro as a unit of length. Then, 
the unit of time will be cRq (or just Ro if c = 1, as in 
the previous equations). The corresponding dimension- 
less variables, normalized by Rq and cRq, will be denoted 
by asterisks. 

So, let us take, for example, r* = 0.01, Rq = 1 and 
study the characteristic shapes of the test-particle orbits 
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FIG. 1: Orbits of the test particles at the specified Schwarz- 
schild radius r* = 0.01 and various de Sitter radii rX (i.e., 
various values of the A-term). 



at various values of r\. Results of numerical integration 
of the equations of motion (f7a)) -(f7c|) for a few values of 
de Sitter radius are shown in Fig. [TJ If r* A = oo (i.e., 
A = 0), the orbit is closed, as should be evidently ex- 
pected. Next, when the values of r A decrease down to 
a few thousand (i.e., the values of A increase), the or- 
bits become slightly spiral; and such an unwinding is ex- 
pressed very well, for example, at r\ — 1000. 

The same orbital radii r* are presented as functions of 
the coordinate time t* (which differs, in fact, only slightly 
from the proper time of the moving particle) in Fig. m 
The curves are wavy because the initial unperturbed or- 
bit was taken to be slightly elliptic. The straight lines 
in this figure represent the standard Hubble motion ex- 
perienced by the test particle (i.e., without any central 
gravitating body) at the same values of r* A . 

It is evident that under certain circumstances (depend- 
ing on the ratio between the above-mentioned character- 
istic parameters) the perturbation of the planetary orbit 
by the dark energy can be appreciable, and the rate of 
increase in the orbital radius can be even comparable to 
the rate of the standard Hubble flow at infinity. In our 
opinion, this points to the potential importance of the 
local cosmological influences on the planetary dynam- 
ics, although a much more careful analysis (with realistic 
planetary parameters and the additional factors affecting 
the planetary dynamics) is still to be done. 
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FIG. 2: Radii of the orbits as functions of time at the specified 
Schwarzschild radius of the central body r'g = 0.01 and various 
de Sitter radii r\ (wavy curves) vs. radii of the test particles 
in the standard Hubble flow (straight lines). 



It should be mentioned also that there is some em- 
pirical evidence in favor of the probable cosmological 
expansion in the Earth-Moon system: This is a well- 
known disagreement between the rates of secular in- 
crease in the lunar semi-major axis measured, on the 
one hand, immediately by the lunar laser ranging [l7j |. 
-Rllr = 3.8 ± 0.1 cm/yr [l8j], and, on the other hand, 
derived indirectly from the data on the Earth's rotation 
deceleration R TOt = 1.6 ±0.2 cm/yr [e.g., the time series 
compiled in [19J . It is quite surprising that these two val- 
ues can be reconciled with each other quite well if, along 
with the commonly-considered tidal interaction between 
the Earth and Moon, one takes into account also a con- 
tribution from the local Hubble expansion; and the rate 
of such an expansion turns out to be in reasonable agree- 
ment with the large-scale cosmological data ■ 

On the other hand, it is commonly believed that the 
modern solar-system ephemerides are able to explain all 
the planetary motions without any additional cosmologi- 
cal influences, apart from the well-known post-Newtonian 
corrections (see, for example, papers [2l|, [22j and refer- 
ences therein). If this is really the case, and no further 



corrections for the A-term are necessary, then the high- 
accuracy planetary observations can be used to impose 
strong constraints on the amount of the dark energy rep- 
resented by the perfectly uniform and constant A-term. 
So, the models with the "dynamic" dark energy (scalar 
field) and/or the nontrivial equation of state may become 
preferable. 

In summary, we presented a rigorous mathematical for- 
mulation of the two-body problem in the A-dominated 
cosmology with an adequate asymptotics at infinity. A 
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set of solutions of the respective equations of motion was 
obtained numerically. It was found that, in some cir- 
cumstances, the secular perturbation of orbits by the A- 
term can be appreciable and even reach the rate of the 
standard Hubble flow. This fact implies that a possible 
presence of the local cosmological influences should be 
taken into account very carefully in the future analyses 
of dynamics of the planetary and other binary systems. 
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